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Abstract. We present results of an analysis of broadband spectra of 14 intermediate polars obtained with the 
RXTE observatory (PCA and HEXTE spectrometers, 3-100 keV). By means of our calculations of the structure 
and the emergent spectrum of the post-shock region of intermediate polars we fitted the observed spectra and 
obtained estimates of the white dwarf masses. The estimated masses are compared with masses obtained by other 
authors. 
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1. Introduction 

Intermediate polars (IP) form a sub-class of magnetic cat- 
aclysmic variables (mCV) (Warner 199.5} . A white dwarf 
(WD) accretes matter from a companion star (typically a 
red dwarf) that fills its Roche lobe. This matter forms an 
accretion disc inside the Roche lobe of the primary which, 
however, is disrupted by the WD magnetic field at some 
distance from the WD surface. As a result the accreting 
matter freely falls on to the WD surface and forms a strong 
shock near its surface. The post-shock matter has a high 
temperature (~ 10-20 keV) and emits X-rays via optical 
thin thermal plasma emission (Lamb & Masters 1979). 

The temperature of the post-shock matter depends on 
the WD mass and therefore the X-ray spectra of IPs can be 
used for WD mass determination (Rothschild et al. 11981(1 . 
Ishida ifTMTl) 



by Aizu l(TrJ75jl . Wu et al. ijigMjl . Woelk & Beuermann 
(1996), Cropper et al. 1|1999[> . see also Frank, King k, 



had estimated the masses of the WDs of IPs 
and polars by using GING A/LAC observations. He fitted 
the X-ray spectra by single temperature bremsstrahlung 
model. Then WD masses were estimated by equaling the 
resulting best fit bresstrahlung temperatures to the max- 
imal shock temperatures. In reality a significant part of 
the post-shock region has smaller temperatures due to ra- 
diative cooling and, as menthioned by the author, WD 
masses obtained in this way are lower limits only. 

In order to estimate the masses of WDs more carefully 
by using X-ray spectra of IPs one needs to calculate the 
temperature and emissivity distribution in the post-shock 
region (PSR). The structure of the PSR was investigated 
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Raine l|2002|l . Wu et al. lfT§9"Hl and Woelk & Beuermann 
(1996} took into account cyclotron cooling, which can be 
important for polars. In these papers it was shown that 
cyclotron cooling is not important for IPs if the surface 
magnetic field is less than rj 10 MG. The model of Wu et 
al. l(T3§Hl was used by Cropper et al. 1(13551 1T555}) for fit- 
ting the spectra of IPs observed by GING A/LAC and by 
Ramsay ( 20001 f° r fitting the spectra of IPs observed by 
RXTE/PCA. Beardmore et al. used the results of 

a simple analytical model of the PSR (see Frank, King & 
Raine 2002) for a comparison with the spectra of V1223 
Sgr observed by GINGA/LAC and ASCA. Ezucka and 
Ishida (1999) estimated the WD masses of nine IPs using 
emission line ratios from ASCA observations. 

The main systematic uncertainty in these works is 
caused by the fact that neither the GINGA/LAC nor the 
RXTE/PCA detectors provide adequate information at 
energies >20-25 keV. Therefore, it is hard to get solid es- 
timates for the masses of the WDs if the mass is > 0.6 M Q , 
because in that case the maximal post-shock temperature 
is > 20 keV. 

In this work we present results of WD mass estimates 
for 14 bright IPs, obtained from broad-band (3-100 keV) 
RXTE spectra. Most of the sources analized have a sta- 
tistically significant signal at energies >10 keV. The ob- 
served spectra have been fitted by PSR model spectra. 
We compare the WD masses resulting from our best fit 
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parameters with the WD masses that were obtained by 
other methods. 

2. The model 

The structure of the stationary post-shock region in plane- 
parallel one-dimensional geometry is described (see e.g. 
Cropper et al. I1999jl by the mass continuity equation 



Tz ipv) 
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Here z is spatial coordinate (see Fig.l), v is the veloc- 
ity of matter, p is the density, T is the temperature, P 
the gas pressure, 7=5/3 the adiabatic index, and p—0.62 
the mean molecular weight of fully ionized matter of solar 
composition. The cooling rate A due to thermal optically 
thin radiation is: 



A 



P 



pm H 



Ajv(T), 



(5) 



where Ajv(T) is the cooling function, here taken for so- 
lar chemical composition as calculated and tabulated by 
Sutherland and Dopita (1993). 
Equation (JTJ has the integral 



pv = a, 



(6) 



where a is the local mass accretion rate per unit area (of 
dimension g cm -2 s _1 ). The local mass accretion rate is 
a free parameter of the model. However, we found that 
the spectrum of emergent radiation is rather insensitive 
to this parameter within reasonable limits (i.e. as long as 
the emission region is optically thin with respect to free- 
free scattering). 

Equations and 10) can be rewritten using JJJJ with 
the substitution z' = zq — z, where zq is the shock coordi- 
nate (see Fig.[Q: 
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dz' 7P — av 



where 



g(z') 



GM wd 
(zo ~ z'? ■ 



(9) 



Equations Q and Q£) were solved by the shooting 
method from z = zq (z' = 0) to z = R wc i {z' = zq — R w d) 
with boundary conditions at z = zq: 



= 0.25V2GAW 
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The WD radius was calculated from the Nauenberg 
l|1972|l WD mass-radius relation: 



R,,„ r , = 7.8 • 10 s cm x 
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In Fig. 2 we present the temperature and the density 
profiles of the post-shock region resulting from our model 
for a particular set of boundary values. For comparison we 
show the profiles obtained by the simple analytical model 
(Frank, King & Raine EDD21 : 



T(z) = To 



P(z) = Po 
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Zq — R 
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This model is based on the assumption of constant pres- 
sure in the post-shock region. In our model the pressure 
grows towards WD surface, therefore the temperature and 
the density in our model are larger. 

As the RXTE/PCA detectors have relatively poor 
spectral resolution and have sensitivity only at energies 
higher than 3 keV we limited ourself to the study of the 
continuum emission. The emergent lines at energies 6-7 
keV in the resulting spectrum were mimicked by broad 
gaussian lines. 

The emergent model spectrum is calculated by sum- 
ming of the local bremsstrahlung spectra: 



Fe = I j(z)dz, 



(16) 



where the local spectra were taken in the following form 
(ZombeckHnnOJ: 



( 8 ) j(z) = 9.52 • 1(T 38 x 
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Fig. 1. . Geometry of a post-shock region model. 



Q. 



t— 







i 










M WD = 0.7 MO 






T/T mlu , 


R WD = 0.01 12 RO 








a = 1.0 g cm -2 s _1 . 




/ 




H s =0.018 R WD 




1 

1 

\ 




H™=0.026 R WD 




\ 








v 












P/Pmax 






, , , 


, i i i i I i i i 



0.005 0.01 
( z - r »d)/ r to 



o.oie 



Fig. 2. . Temperature and. density profiles for one of the 
PSR models. Dashed lines renote simplistic analytic solu- 
tion of the structure of post-shock region (Frank, King & 
RaineHHH 



3. Observations and data reductions 

In order to compare the results of our model calculations 
with observations we have chosen the data of the RXTE 
observatory that cover a relatively broad energy band (3- 
250 keV), where IPs release most of their emission. We 



limited ourselves only to relatively bright IPs, as we are 
interested in the hard X-ray (>20 keV) part of the spectra 
of IPs. The list of IPs consists of V1223 Sgr, FO Aqr, 
EX Hya, AO Psc, TV Col, GK Per, V709 Cas, PQ Gem, 
V2400 Oph, BG CMi, V405 Aur, V1062 Tau, DO Dra, 
TX Col. 

The data of the RXTE/PCA and the RXTE/HEXTE 
were reduced by means of standard tasks of the 
LHEASOFT /FTOOLS 5.2 package. The background sub- 
traction of the PCA detectors was done with the CM17_240 
model. The subtraction of the background of the HEXTE 
detectors is especially important in our case which con- 
centrates on the hard X-ray part of the spectrum of IPs. 
In order to monitor the background of the detectors the 
HEXTE is rocking by ±1.5° with respect to the target. 
We analised the HEXTE spectra of IPs only if the back- 
ground spectra obtained at the +1.5° and the —1.5° offset 
positions were statistically compatible, in other words if 
there were no contaminating sources in the vicinity of our 
targets. 

4. Results 

The spectra of the 14 IPs, extracted from the RXTE data 

are presented in Figs. El El 

For fitting these spectra we used results of our model 
calculations described in the previous section with one or 
two (only in the case of very bright IPs V1223 Sgr and GK 
Per) partial covering components. In contrast to the model 
used by Cropper et al. Ijl998|l our model lacks the reflected 
component. However, as has been shown by Cropper et 
al. Ijl998|l taking into account the reflection changes the 
best fit parameters and thus the WD mass only slightly. 
In addition to that, the usage of spectral information at 
energies higher than ~20 keV helps us to further reduce 
the influence of the reflected component on the shape of 
the spectral cutoff. 

The best fit parameters resulting from our model are 
presented in Table Using the observed and absorption- 
corrected fluxes of IPs we esimated their luminosities and 
mass accretion rates. The best fit model is shown in Figs. 
El - El by a solid line. 

We compare the resulting estimates of WD masses 
with those measured by other methods. In the catalog 
of Cataclysmic Variables (Ritter & Kolb 2003) there are 
only 5 intermediate polars that have a spectroscopically 
measured mass. On Fig.EJwe present a comparison of our 
mass estimates with those given in the catalog. It is seen 
that we have a relatively good agreement, however the er- 
ror bars of all these estimates are quite large. The only 
system that likely has a higher mass than we obtained 
from our X-ray spectral modelling is GK Per (during out- 
burst). One possible explanation of this fact could be the 
WD's fast rotation (P sp = 351 s). It is well known that 
the radius of the magnetosphere in systems with an accre- 
tion disc depends on the mass accretion rate. If the WD 
in GK Per is rotating close to its equilibrium period in 
its quiescent state (i.e. if the radius of its magnetosphere 
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Table 1. Best fit model parameters and inferred physical system parameters. Here Cf the partial covering coefficient, 
F3-100 the observed X-ray flux in the band 3 - 100 keV, Fo. 1-100 the unabsorbed X-ray flux in the band 0.1 - 100 
keV, both in 10 -11 erg cm -2 s _1 units. Distances for V405 Aur, V709 Cas, PQ Gem and V2400 Oph are assumed. 



Name 


M wd , 


C F 


Nh, 


F3-100 


Fo. 1-100 


M 


La; 


d 




M Q 




10 22 cm~ 2 






10 16 g s" 1 


10 33 erg s _1 


pc 


V1223 Sgr 


0.95 ±0.05 


1.0° 


3.7 ±1 


4.2 


7.3 


12.6 


24.3 


527 1 






0.57 ±0.05 


28 ±3 












FO Aqr 


0.6 ±0.05 


0.9 ±0.1 


36 ±4 


1.5 


4.9 


10.2 


9.4 


400 2 


EX Hya 


0.5 ±0.05 


1.0 


0.6 ± 1.0 


1.5 


3.4 


0.25 


0.17 


64.5 3 


AO Psc 


0.65 ±0.05 


0.95 ±0.1 


10 ±5 


1.3 


2.5 


4.7 


4.9 


420 4 


TVCol 


0.84 


0.81 


21 


1.7 


3.5 


3.4 


5.7 


368 5 


GK Per 


0.59 ±0.05 


1.0 a 


20 ±2 


12.5 


88.6 


81.5 


156 


340 6 






0.75 ±0.05 


125 ± 10 












V709 Cas 


0.9 ±0.1 


0.5 ±0.1 


29 ±5 


2.3 


3.7 


5.76 


11.1 


(500) 


PQ Gem 


0.65 ±0.2 


1.0 


5±1 


0.80 


1.3 


3.7 


3.9 


(500) 


V2400 Oph 


0.59 ±0.05 


0.4 ±0.1 


29 ± 10 


2.43 


5.6 


18.7 


16.75 


(500) 


BG CMi 


0.85 ±0.12 


0.70 ±0.05 


25 ±6 


0.77 


1.54 


5.3 


9.0 


700 7 


V405 Aur 


0.9 ±0.1 


0.6 ±0.1 


39 ±3 


0.9 


1.7 


2.65 


5.09 


(500) 


V1062 Tau 


1.0 ±0.2 


0.6 ±0.1 


21 ±7 


1.1 


1.8 


2.2 


5.4 


500 s 


DO Dra 


0.75 ± 0.05 


0.4 ±0.1 


20 ±5 


7.8 


14.4 


3.1 


4.1 


155 9 


TX Col 


0.7 ±0.3 


< 0.4 


<60 


0.32 


0.58 


1.8 


2.1 


550 10 



a - for systems V1223 Sgr and GK Per which have very high mass accretion rates we used two absorbers. 

References to distance estimates: 1. Beuermann et al. (2004), 2. McHardy et al. 1987 3. Beuermann et al. 2003 4. Hellier, 
Cropper & Mason UM!^ . McArtur et al. l2l7fTT1 6. Warner IM 7. Berriman ITiWl 8. Szkody & Silber ITiMl 9. Mateo, Szkody 
& Garnavich [TMn 10. Buckley & Tuohv H9l&l 



Table 2. Comparision of mass estimates obtained in this 
work with those obtained with the usage of other X-ray 
instruments 



Name 


PCA±HEXTE 


RXTE/PCA 


Ginga 


ASCA 




M e 


M Q 


M e 


Me 


V1223 Sgr 


0.9 


1.1 




1.28 


FO Aqr 


0.6 


0.88 


0.92 


1.05 


EX Hya 


0.5 


0.45 


0.46 


0.48 


AO Psc 


0.65 


0.60 


0.56 


0.40 


TV Col 


0.84 


0.96 


1.30 


0.51 


GK Per 


0.59 






0.52 


V709 Cas 


0.9 


1.08 






PQ Gem 


0.65 








V2400 Oph 


0.59 


0.71 




0.68 


BG CMi 


0.85 


1.20 


1.09 




V405 Aur 


0.9 


1.10 




>0.54 


V1062 Tau 


1.0 


0.90 






DO Dra 


0.75 








TX Col 


0.7 


0.73 


0.48 


0.66 



Therefore the mass of the WD could be underestimated 
by as much as 10-20%. 

Note that our mass estimates are in better agreement 
with those obtained from optical methods than those ob- 
tained by Cropper et al. (1999) from their analysis of the 
GINGA/LAC X-ray spectra (-2-20 keV) (see Table©. 
The masses obtained by Ramsay [2000 from his analysis 
of the RXTE/PCA (-2-20 keV) spectra are closer to our 
estimates, but also slightly higher. A possible reason for 
this difference could be the lack of the hard X-ray part in 
the observed spectra. 

It is interesting to mention that the parameters in our 
Table show a correlation of the mass accretion rate with 
the strength of photoabsorption, that, in principle, is an- 
ticipated for IPs because systems with higher mass ac- 
cretion rates have higher Nh values. A similar effect was 
observed in the case of GK Per (quiescent and outburst 
spectra) by Ishida et al. I)1992|l . 



is close to the corotation radius) then it means that dur- 
ing an outburst (increase of the mass accretion rate by a 
factor of 10-20, see Ishida et al. H992f) the size of the mag- 
netosphere should shrink by a factor of — 2. Therefore the 
radius of the magnetosphere in GK Per during an out- 
burst can be as small as 2 — 3 x 10 9 cm, and thus, only a 
factor of 5-8 larger than the radius of the WD. In this case 
the free fall velocity gained by the accreting matter near 
the surface of the WD will be somewhat smaller than in 
the case of falling from infinity (our model assumption). 



5. Conclusions 

Estimating the masses of the WDs in cataclysmic variables 
is a relatively complicated problem. The standard method 
based on radial velocity measurements often gives unreli- 
able results because of our poor knowledge of the incli- 
nation angle of the binary system. Besides, spectral lines 
from accretion discs are very broad in comparison with 
Dopplcr shifts caused by orbital motion and their cmissiv- 
ity distributions can be non-axisymmetric (in the frame of 
reference of the WD). Therefore, additional assumptions, 
such as a relation between the mass of the secondary star 
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Fig. 6. Comparison of the masses of the white dwarfs 
as estimated from the RXTE broadband spectra (this 
work) with estimates obtained by other methods. (1)- 
Beuermann et al. (2003), (2) - Crampton et al. (1986), (3) 
- Reinsch (1994), (4) - Haswell et al. (1997), (5) - Mateo 
et al. (1991) , (6) - Friend et al. (1990), (7) - Hellier (1993) 
, (8) - Penning (1985) 



and the orbital period (based on a mass-radius relation 
for main sequence stars, see Warner H.995[l are often used. 

Magnetic cataclysmic variables provide us with an ad- 
ditional way for estimating the masses of the WDs. The 
matter accreted onto the WD surface is heated in a stand- 
ing shock and the maximal temperature of the plasma af- 
ter the shock front depends practically on WD mass only. 
Therefore it is possible to estimate WD masses in polars 
and intermediate polars by fitting their X-ray spectra with 
model spectra of the radiation emerging from the post- 
shock region. If both methods of mass determinations are 
available it is very interesting to compare the resulting 
mass estimates. The temperatures of the PSR are rela- 
tively high (10 - 30 keV). Therefore the study of WD 
spectra in standard X-ray band suffers from a lack of in- 
formation at energies of a few x kT, that would strongly 
help to constrain the WD mass estimate. 

In this work we present the WD masses of 14 
IPs estimated from fitting the broad band (3-100 keV) 
RXTE/PCA+HEXTE X-ray spectra with model spectra 
of the PSR. We did not take into account cyclotron cooling 
of the PSR and therefore did not consider polars, where it 
can be important. Our estimates of WD masses are smaller 
than those obtained from GING A/LAC (Cropper et al. 
HMty and RXTE/PCA CRamsav lMTUjl data. We also find 



satisfactory agreement between our WD mass estimates 
and those obtained by optical spectroscopic methods. We 
conclude that it is necessary to include the hard X-ray 
spectral band (E > 20 keV) for obtaining reliable WD 
mass estimates in mCVs. 
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Fig. 3. Spectra of intermediate polars obtained with the RXTE observatory. Crosses denote the PCA data, open 
circles - HEXTE data. Solid lines show best fit model with parameters presented in Table 1 
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Fig. 5. The same as Fig. |3 



